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Abstract 

We describe the behavior as n — > oo of the Laplace transforms of P n , where P a 
fixed complex polynomial. As a consequence we obtain a new elementary proof of an 
result of Gillis-Ismail-Offer |2j in the combinatorial theory of derangements. 

1 Statement of the main results 

The generalized derangement problem in combinatorics can be formulated as follows. Sup- 
pose X is a finite set and ~ is an equivalence relation on X. For each x £ X we denote by 
x the equivalence class of x. X^ will denote the set of equivalence classes. The counting 
function of ~ is the function 

v = v^, '. X — > Z, u(x) = \x\. 
A ~- derangement of £ is a permutation <p : X — > X such that 

x x, Vx £ X. 

We denote by ~) the number of ~-derangements. The ratio 



\X\\ 



is the probability that a randomly chosen permutation of X is a derangement. 

In S. Even and J. Gillis have described a beautiful relationship between these numbers 
and the Laguerre polynomials 

k=0 v 7 

For example 

Lq(x) = 1, L 1 {x) = l-x, L 2 {x) = ^(x 2 -4x + 2). 
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We set 

L~ := X\(-iy {c) v{c)\L v[c) (x). 

cex 

Observe that the leading coefficient of is 1. We have the following result. 
Theorem 1.1 (Even-Gillis). 

/•oo 

W(X,~) = / e~ x L^(x)dx. (1.1) 
Jo 

For a very elegant short proof we refer to [3]. 

Given (X, ~) as above and n a positive integer we define (X n ,~ n ) to be the disjoint 
union of n-copies of X 

n 

X n ={JXx{k} 

k=l 

equipped with the equivalence relation 

~ n (y, k) <^=>- j = k, x ~ y. 

We deduce 

1 f°° 

p{Xn ~ n) = J^\x\y. L e ~ x ( L ~^ ndx ^ 

For example, consider the "marriage relation" 

(C,~), C = {±1}, -1-1. 

In this case C consists of a single element and the counting function is the number v = 2. 
Then (C n j ^72) Cciii be interpreted <is ci §roup of ?r iriflxriecl couples. If we set 

5 n :=p(C n ,~ n ) 

then we can give the following amusing interpretation for S n . 

Couples mixing problem. At a party attended by n couples, the guests were asked to 
put their names in a hat and then to select at random one name from that pile. Then the 
probability that nobody will select his/her name or his/her spouse's name is equal to 5 n . 

Using (|1.2|) we deduce 

We can ask about the asymptotic behavior of the probabilities p(X n ,~ n ) as n — > 00. In 
0, Gillis- Ismail-Offer describe the first two terms of an asymptotic expansion in powers of 
71 . To formulate it let us introduce "momenta" 

v r = ^v{c) T . 

cex 
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Theorem 1.2 (Gillis-Ismail-Offer). 

P (X n , ~ n ) = eM-~) ( 1 " ^-f- V K -' + 0(n- 2 ) as n ^ oo. (1.4) 
For example we deduce from the above that 

S n = e~ 2 (l - -n- 1 + 0(n- 2 ), n -» oo (1.5) 
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The proof in [2] of the asymptotic expansion (|1.4|) is based on the saddle point technique 
applied to the integrals in the RHS of (|1,2|) and special properties of the Laguerre polyno- 
mials. 

In this paper we will investigate the large n asymptotics of Laplace transforms 

dn+l poo 

3>i(9> z) = tj^t / e- zt Q(t) n dt, fttz > 0, (1.6) 
(dn)\ J 

where Q(t) is a degree d complex polynomial with leading coefficient 1. If we denote by 
£j[f(t),z] the Laplace transform of f(t) 



/•oo 

L[f(t),z]= / e~ zt f{t)dt then 9 n (Q,z) 
Jo 



H[Q(t) n ,z] 



L[t dn ,z\ 

The estimate (|1.4j) will follow from our results by setting 

z = l, Q = L_ 

To formulate the main result we first write Q as a product 

d 

i=i 

We set 

1 d 

r = {r u ---,r d )£C d , /x s = n s {r) = - ^Vf 

i=l 

Theorem 1.3 (Existence theorem). For every 9^e z > we have an asymptotic expan- 
sion as n — > oo 

oo 

3 n {Q,z) = {Y,A k {z)n- k ). (1.7) 

k=0 

Above, the term Ak(z) is a holomorphic function on C whose coefficients are universal 
elements in the ring of polynomials C(d)\pi, fi2, • • • , fJ>k}> where C(d) denotes the field of 
rational functions in the variable d = degQ. 

We can say a bit more about the coefficients Ak{z). 

Theorem 1.4 (Structure Theorem). For any k and any degree d we have 

A k (z) = e^ z B k (z), 

where B k £ C(d)\p\, ■ ■ ■ , fj>k][z] is a universal polynomial in z with coefficients in 
C(d)[fj,i, ■ ■ ■ 
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Here is a brief description of the contents. We give the proof of the existence theorem 
in the next section, while in the third section we compute the terms in some cases. For 
example we have 



A (z) = e^\ A x {z) = ^ z (pj - n 2 )z\ 



and we can refine ()1.5|) to 

/ 1 23 

5 n = e- 2 (l--n- 1 -—n~ 2 + 0(n' 2 ), n -> oo. (1.8) 

These computations will lead to a proof of the structure theorem. 

For the reader's convenience we include a list of symbols we will use throughout the 
paper. The symbol [n] denotes the set {1, 2, • • • , n}. A d- dimensional (multi)index will be 
a vector a G 2>o- For every vector x G C d and any <i-dimensional index a we define 

x a = x" 1 ■ ■ ■ x^ d , \a\ = a\ + ■ ■ ■ + ay, S(x) = x\ H + x^. 

If n = \a\ then we define the multinomial coefficient 

f n\ n! 



x aj Y[ i=1 aiV 
These numbers appear in the multinomial formula 



— H\ 

XT. 



\a\=n 

2 Proof of the existence theorem 

The key to our approach is the following elementary result. 

Lemma 2.1. If P(x) = p m t m + • • ■ pit + po is a degree m with complex coefficients then for 
every JHc z > 1 we have 

1 ' J J[) a+b=m \aJ 



Proof 



'in+1 poo „m+l m />oo 'm+1 m „i 

7 e~*P{t)dt= Z —Y,Pa ^~E^ir= E 7& A 



m! 



m ) 6! 

a / 



a=0 a=0 a+b=m 

□ 

Denote by Q(n,a) the coefficient of t a in Q(t) n . From ()2.1() we deduce 

w- E <"> 
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Using the equality 

d / n 

Qn =n ( e cy 

i=l \j+k=n 



, fc 



(t+r-i) n 

we deduce that if a + 6 = (in then 



Q(n,a)=E (n(r))^ (2 - 3) 

|5|=6 \i=l V 3/ J 



For |a| = b we set 

Li W © 

so that 



J n (Q,z) = E P ^(«)M«))-^ (2-4) 

a+b=dn \\a\=b / 



Observe that we have 



p ... nt,p-i)-a-«g t ) i/i'. 



o 



The coefficients p,(a) define the multinomial probability distribution Pj, on t ne se t °f rnul- 
tiindices 

A t = {ae (Z> ) fe ; |a| =&} 

For every random variable £ on A& we denote by Pf,(C) its expectation with respect to the 
probability distribution P&. For each n we have a random variable Cn,b on A& defined by 

, f ^ ntia-^)---(i-^) 

Cn,6(«) = =T^T- j" M«)- 

Form (|2.4|) and (|2.5() we deduce 



? n (Q,z)= E b((n,b)^- (2-6) 

a+b=dn 

To find the asymptotic expansion for 3~ n we will find asymptotic expansions in powers of 
nT 1 for the expectations Pfe(Cn,&) and them add them up using (|2.6[) . 
For every nonnegative integer a we define a polynomial 



W a (z) = { 



if a = 0,1 



i nSa-^) if «>i 
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For a d-dimensional multi-index a we set 

d 



W s (x)=l[W CH (x). 



i=l 



We can now rewrite (|2.5[) as 



W^ 1 



P n>h (a) = P b {a)-^t L 



We set 

Rb(a,x) = Ws(x), K b (a,x) = — R b (a, x)p b (a). 
We regard the correspondences 

a i— ► R b (d, x), K b (d,x) 

as a random variables R b (x) and K b {x) on A& valued in the field of rational functions. We 
deduce 

Cn,b = K^n" 1 ). 

Observe 

E b (x) = E b {K b {x)) = ^—E b {R b (x)) 

From the fundamental theorem of symmetric polynomials we deduce that the expectations 
E b (R b (x)) are universal polynomials 

E b {R b (x)) £ C[m, ■ ■ ■ ,/J, b ][x], deg x E b (R b (x)) <b-d, 

whose coefficients have degree b in the variables /Xj, deg^j = i. We deduce that E b (x) has 
a Taylor series expansion 

E b (x) = E b (m)x m 

m>0 

such that E b {m) £ C(d)[pi, ■ ■ ■ ,p b ]. The rational function x — > K b (a,x) has a Taylor 
expansion at x = convergent for \x\ < r^r so the above series converges for \x\ < t^t. We 
would like to estimate the size of the coefficients E b (m). The tricky part is that the radius 
of convergence of E b {x) goes to zero as b — > oo. 



Lemma 2.2. Set 



R = max I r j I . 

Ki<(i 



There exists a constant C which depends only on R and d such that for every b > and 
every 1 < X b < rrj we have the inequality 

h \ m c b 

\ b d) (6-2)!(l-A 6 ^i)' 
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Proof Note first that 

\p b {a)\ < R b , V|a| = b. 

For b = 0, 1 we deduce form the definition of the polynomials W a that E b {x) = 1. Fix m 
and b > 1. Using the Cauchy residue formula we deduce 

E b ( m ) = 1 ^ i J^— l E b { X )d X , H = \ b --. 

Hence 

1 b m R b 

\E b (m)\ < — sup 1-^6(^)1 < 7T ,x m — : , w 7 ,,x, ■ max E b (R b (x)) . 

fr m \x\=h {\ h d) m min\ x \ =h \W b {x / d)\ \ x \=h 

Next observe that 

6-1 

W 6 (a;/d) = (b - 1)1 Y[(r- x/d)), h/d < l/k,Vk <b-l, 



"k 
k=i 



from which we conclude 



min\W b ( X )\=W b (h)Jli(l-^)='li(b-kX b )> {b 2)K1 V ^ 



1 1 k=l k=l 

To estimate E b {R b {x)) from above observe that for every 1 < k < (b — 1) and |x| = h we 
have 

, , , fcAftd 
|1 - jfea?| < 1 + k\x\ = 1 + — < 1 + d 

This shows that for every \6t\ = b and |x| = h we have 

|i? 6 (a,x)| < (l + d) b . 

The lemma follows by assembling all the facts established above. 

□ 

Define the formal power series 

z b 

A m {z) :=J2 E b(m)-eC{z}. 

b>0 

The estimate (|2.7j) shows that this series converges for all z. 

For every formal power series / = X^fc>o a kT k and every nonnegative integer £ we denote 
by 4(/) its i-th jet 

4(/) = E Q * Tfc - 

k=0 

For x = n _1 we have 

6<d/a: ' 6<d/x m>0 
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E ( E ^M^)* m = E ^'(^m*)) 

m>0 6<d/x ' m>0 



x m . 



Consider the formal power series in x with coefficients in the ring C{z} of convergent power 
series in z 

9oo(V) = E A m{z)x m € C{z}[[x]]. 

m>0 

We will prove that for every t > and every z G C we have 

|J„(z) - 49oo(V)| = 0(n~* _1 ), asn^oo. (2.8) 
To prove this it is convenient to introduce the "rectangles" 

D u , v = {(b,m) <E (Z> ) 2 ; b < u, m < u.j 
In this notation we have (x = n _1 ) 

•J n {z)= E E b(m)x m ^, 4?co(z)= E ^M^IT- 



6!' ^ -uv-v- 6| 



Then 



- j^oo(^) = EE^( m ^ m b + E E E *( 



u 1 ^v^"" lmj ¥. |a 

b<dn m>£ m<l b>dn 



Si(n) S 2 (n) 

We estimate each sum separately. Using 1)2.7(1 with a A;, > 1 to be specified later we deduce 

> LB 6 (m)x m < in-> (^— ;) • 

1 " (6-2)!(l-A 6 ^i) ^/A^/ 

The inequality b < dn can be translated into *j < 1 so that the above series is convergent 
for b < dn whenever A^ > 1 so that 

|v^ ln/ . C b b h - X / bx \t+i 1 



(6-2)!(l-W)VA 6 d; i - to- 



When b < dn we have 

fee 1 
Aftd A;, 



If we choose 



& \l/2 



6-1 



we deduce 

6-1, /6-l\i/2 1 



1-A6— ) = l-( — ] => - , b -i <b 

1-A 6 — 
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and, since ^ < 

1 1 
7— < r < 2b. 

1 ox i _l_ 

Using the inequalities 

-J , VA; > 

we conclude that for 6 < dn we have 

^ |£ 6 (m)x m | < C^V^. 

Since the series $^&>o^I^ +2 fr conver g es we conclude that 

S x {n) = 0(x m ). 

To estimate the second sum we choose = 1 in (|2.7j) and we deduce 

E h (m) < Cl 

Hence 

\X-F( ^\^ {CM) b b 2 ,.. r ■ (|C 3 ^|) b - 2 
|^^(m) ¥ |< <(2C 3 N) L (5-2)! 

Using Stirling's formula we deduce that for fixed z we have 

b>dn y 1 

Hence 

|5 2 (n)| < C 4 (z){£ + l)n- e -\ 
This completes the proof of (|2.8j) and of Theorem 11.31 



3 Additional structural results 
§3.1 The case d = 1. Hence 

Q(i) = (t + Ml)- 

We have 



/■°° /*°° 77I 

/ e-**(i + m) n dt = e^ z / e - rf t n dt = e^—-^ 
Jo Jo z 



Hence in this case 
and we deduce 

A (z) = e wz , A k (z) = 0, Vk > 1. 



□ 
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§3.2 The case d = 2. This is a bit more complicated. We assume first that /ii = so 
that 

Q(t) = t 2 - a 2 



Then 



and we deduce 



Q(n, a) = < 



f (-l)V 2 ("- fe )(£) if a = 2fc 
if a is odd 



rr , v y (~l) fc L" b ) (^) 2fe y n\{2n - 26)! (-1) W fe 
^ h ( 2 n- 26 ) « fe(»-fr)K2n)! 6! 

_A n(n-l)---(n-b+l) (-l) b (az) 2& 
~ ^ 2n(2n - 1) • • • (2n - 26 + 1) 6! 

A J_ _ 6 (l-l/n)---(l-(ft-l)/n) (-l)W 
Z^2 2b (1 - l/(2n) ■ • • (1 - (26 - l)/(2n) 6! 

i _ i-i i M2 , 1-2 (i - iA0 Ml , 

2 l-l/(2n) 2! ^2 4 (1 - l/(2n))(l - 2/(2n))(l - 3/(2n)) 4! 

To obtain A^{z) we need to collect the powers n~ k . The above description shows that the 
coefficients of the monomials z 2b contain only powers n~ fe , k > 6. We conclude that A^{z) 
is a polynomial and 

deg z A k (z) < 2k. 
Let us compute the first few of these polynomials. We have 

f „ ( ,, = i-i„->(i + i„- + ...)M! + ^(i + ...)^ + ... 

We deduce 

A (z) = l, A l (z) = -\(az) 2 1 A 2 {z) = - l -{oz) 2 + ^{az)\ 
If /xi 7^ so that 

Q(t) = (t + ri)(t + r 2 ), ri+r 2 = 2/ii, 
then we make the change in variables t = s — m so that 

Q(t) = P( S ) = s 2 - r 2 , a 2 = (n - W) 2 = (ri " 4 r2)2 . 

Now observe that 

4/j 2 + (n - r 2 ) 2 = (n + r 2 ) 2 + (n - r 2 ) 2 = 2(r 2 + r£) = 4,u 2 

so that 

2 9 

a =/x 2 -Mi- 
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Then 

2n+l roo „2n+l roo 

3" n (Q, z) = -— \ e- zt Q(t) n = — — / e-< s -^P(s) n ds = e^ 12 J n (P, z). 
We deduce 

A (z)=e^, Al (z) = - e —(azf, A 2 (z) = (- -(a zf + -^{a zf) . (3.1) 

For the couples mixing problem we have 

Q(t) = t 2 - At + 2 

so that 

Mi = ~\ = "2, a 2 = \(n- r 2 f = i ((n + r 2 ) 2 - 4 n r 2 ) = 1(16 - 8) = 2. 
and we deduce 

5 n = J n (Q, z = l) = e~ 2 (l - V 1 - ^n~ 2 + 0(n^ 3 )) . (3.2) 

§3.3 The general case. Let us determine the coefficients Aq(z) and A\{z) for general 
degree d. We use the definition 

A k {z)=Y J E b {k)^- 

b>0 

For \a\ = b 

d ai — 1 

w s {x) = w b , a ( x ) = n( n (! - 
i=i j=i 

d / a;— 1 \ 1 rf 

II i (E •/)•'• • ••• =i-2(E a *(°«- 1 )) a; + "- • 

i=i V j=i / i=i 

w 6 (x/d) = na+i^+---) = i+^=^^+--- • 

fc=l 

Next compute the expectation of R b (x) 

1 d 

E b (R b (x)) = E b { Pb ) - ^(Yl " !)^)^ + 

i=i 

The multinomial formula implies 

£ 6 (p 6 ) = Mi- 

Next ^ 

i=l |a|=6 V 7 i=l 
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Now consider the partial differential operator 

d d 2 
i=i * 

Observe that the monomials r are eigenvectors of y 

d 



i=l 



We deduce 



Hence 



and we deduce 



E h {R b {x)=iA-^{'?iii)x + 



E b (x) = U-^ b 1 )x + ...)(l + b -^lx + ---) 



We deduce A (z) = e^ z 

2 00 b -1 2 2 1 

6=2 

We can simplify the answer some more. 



Al ( z ) = ^y - iov** = ^ e w* - iov** 

v ; (6-2)! 2 2d 2 



i=l 

We conclude that 



d (6-2)! d 

b>2 y ' 

Hence 

A (z) = e^, A l {z) = —^\- l x 2 )z 2 
For d = 2 we recover part of the formulae . 
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§3.4 The proof of Theorem 11.41 Clearly we can assume d > 1. We imitate the 
strategy used in the case d = 2. Thus, after the change in variables t — > t — fj,i we can 
assume that /zi = so that Q(t) has the special form 

Q(t) =t d + a d _ 2 t d - 2 + --- + a . 

Set 

T( „, 6):= 5(!^) 

Idn-bJ 

This is a power series in x = n" 1 , 

T(n,b)=T b (x)\ x=n -i, T b (x) =Y,T b {k)x k . 

k>0 

We have 

A k [ Z ) = Y J n{k) z ^ 

and we need to prove that Aj. is a polynomial for every k. We denote by £(b) the order of 
of the first non-zero coefficient of T b (x), 

i(b) = min{£; > 0; T b (k) / 0}. 

To prove the desired conclusion it suffices to show that 

lim i{b) = oo. (3.4) 

For every multiindex = (0 d , d -2, ■ • • ,0i,Po) we set 

L($ = d&+(d-2)fl,_ 2 + ... + /?i. 
We set a = (1, a^_2, • • • , oi, oo) G We have 

nn,6) = -i-. £ (3.5) 

Un-W |/3|=n,L(/3)=dn-6 

Now observe that we have 

d\0 - L0) = 2f3 d _ 2 + 3/3 d _ 3 + • • • (d - l)Pt + dp = d\$\ - L0) = b. 
In particular we deduce 

0i < — — <-, V0 < j <d-2 (3.6) 
' 3 ~ d-j ~ 2 ~ v; 

and 

2/? d + 6 = 20 d = 20 d _ 2 + 3/3 d _3 + • • • (d - 1)A + dft 
> 2/3 d + 2/? d _ 2 + • • • + 2f3 1 + 2/3 = 2n 

so that 

n-0 d < b -. (3.7) 
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These simple observations have several important consequences. 

First, observe that they imply that there exists an integer N(b) which depends only 
b and d, such that for any n > the number of multi-indices (3 satisfying \/3\ = n and 
L(f3) = dn — b is < N(b). Thus the sum ()3.5|) has fewer than N(b) terms. 

Next, if we set \a\ := max < :( .<rf_2 then, we deduce 

\a?\ < \ a f a+ -+^ < |a|^ = C 5 (b). 

Finally, using the identity 

/n\_/n\ (n- f3 d \fn- (3 d - f3 d ^ 2 \ 
W"WV^- 2 A Pd-3 )"' 

the inequalities (|3.7j) and (™) < 2 m , Vm > k we deduce 

®<(£.^^( » )<W»w« v„», 



Hence 



|/3|=n,L(/3)=dn-6 
On the other hand 



£ I (j?K| ^ ^(^)C 5 (6)C 6 (6)n^J + i = cv^L^J+i. 



< C s (b)n 



I dn \ 
\dn-b) 

so that 

\T(n,b)\ = inin- 1 )] < C 9 (&)J 6 /2j+i-f> < C 9 {b)n 1 - b / 2 . 

This shows 

T b (k) =0, Vk < b/2 - 1 

so that 

1(b) > b/2 - 1 -> oo as 6 oo. 

Remark 3.1. We can say a bit more about the structure of the polynomials 

B k (m, ■■■ ,Hd,z) € Rd = C[/ii, • • • /c > 0. 

If we regard B as a polynomial in ri,- •■ , we see that it vanishes precisely when ri 
• • • = r d . Note that 

ri = • • • = r d = r «=► Q(t) = (* + 

On the other hand 

■ k = 1 hh {rtt) = d^T^~t = d2^7T^r~dW)' 

k i=l k>0 i=l i=l 



□ 
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If Q(s) = (s + r) d we deduce 

s Q'(s) s _ 1 
d Q(s) ~ s + r ~ l-rt 

This implies that 

n = ■■■ = r d n\ = fi % j, VI < i, j < k 

The ideal / in Rd generated by the binomials fj,{ — fij is prime since Rd/I = C[//i, z]. Using 
the Hilbert Nullstelensatz we deduce that must belong to this ideal so that we can write 

Bkim,-- ■ ,f*d,z) = A 2 k{^,z)(nl - fi 2 ) H h A dk {n,z)(i4 - fj, d ). 

□ 

References 

[1] S. Even, J. Gillis: Derangements and Laguerre polynomials, Math. Proc. Camb. 
Phil. Soc, 79(1976), 135-143. 

[2] J. Gillis, M.E.H. Ismail, T. Offer: An asymptotic problem in derangement theory, 
Siam J. Math. Anal., 21(1990), 262-269. 

[3] D.M. Jackson: Laguerre polynomials and derangements, Math. Proc. Camb. 
Phil. Soc, 80(1976), 213-214. 



fc>0 



jij = n[, VI < j < d. 
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